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Abstract 

We extract the heavy quark masses and non-perturbative parameters from the Delphi prelim- 
inary measurements of the first three moments of the charged lepton energy and hadronic mass 
distributions in semileptonic B decays, using a multi-parameter fit. We adopt two formalisms, one 
of which does not rely on a l/m c expansion and makes use of running quark masses. The data 
are consistent and the level of accuracy of the experimental inputs largely determines the present 
sensitivity. The results allow to improve on the uncertainty in the extraction of \V c b\. 
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1 Introduction 



The Operator Product Expansion (OPE) represents a foundation for extracting the \V ub \ and |V cb | 
elements of the CKM mixing matrix from inclusive semileptonic (s.l.) B decays. In this framework, 
the decay width is expressed in terms of quark masses, and non-perturbative effects are described 
by expectation values of heavy quark operators, some of which are presently poorly known. The 
experimental accuracy already achieved, and that expected from the large data sets recorded by the 
5-factories, makes the ensuing theory uncertainty a serious limitation. Extracting the heavy quark 
masses and the non-perturbative parameters, arising from the 1/mjj and 1/rajj corrections, directly 
from the data has therefore become a key issue. There have already been \V C ^\ determinations from 
the first moment of distributions in s.l. and b — > X s ^ decays, and the 1/mjj corrections, estimated 
from parameter ranges, have been found to represent an important source of uncertainty [1]. These 
ranges, based on dimensional arguments, are affected by a considerable degree of arbitrariness. 

In order to circumvent these problems, we introduce in this Letter a multi-parameter fit to 
determine the relevant l/mf and l/m\ parameters, together with the heavy quark masses, from the 
first three moments of the leptonic energy and hadronic mass spectra in s.l. B decays. Results are 
based on preliminary data obtained by the Delphi Collaboration. Moments are measured without 
cuts on the lepton energy in the B rest frame. We consider two formalisms, one of which is new 
and relies on fewer theoretical assumptions. The use of higher moments guarantees a sensitivity 
to these parameters and the simultaneous use of the hadronic and leptonic spectra ensures that a 
larger number of parameters can be kept free in the fit. We discuss the results both in terms of 
the extraction of the parameters and the implications for \V c b\, and as a consistency check of the 
underlying theoretical assumptions. 

2 Extracting Non-Perturbative Parameters 

The moments of the hadronic and leptonic spectra in s.l. B decays have recently been measured 
by several experiments [1-5]. We consider here moments of the charged lepton energy distribution 
defined as 

M 1 (E t ) = ± JdE e E e ^ M n (E e ) = i J dE e (E e - M x {E t )) n J^- (n > 1), (1) 
and moments of the distribution of Mx, the invariant hadronic mass, 
M 1 (M X ) = ± \dM\{Ml-Ml)^ M n {M x )= l - f dM x {M x -{M x )T^ (n > 1), (2) 

X X 

where M D = 1.973 GeV is the spin averaged D meson mass and no cut on the charged lepton energy 
is assumed. 

The theoretical framework to interpret these data has long been known and it is based on the 
OPE. Different implementations exist, depending on the way the quark masses are treated. For 
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instance, the mj and m c masses can be taken as independent parameters or subject to a constraint 
on mb — m c , imposed from the measured fiW an d £)(*) meson masses. The second choice introduces 
a l/m c expansion. Another option concerns the normalization scheme used for quark masses and 
nonperturbative parameters. One approach is to use short-distance masses, such as the low-scale 
running masses. Alternatively, the pole mass scheme can be used. 

The OPE expresses lepton moments through quark masses as a double expansion in a s and l/m&: 



\I,,(E,)= 



Mr) + <in(r) - + b n (r) 4 + Ur) 4 + Ur) 4 + 4f + •■ 
7r m 6 ml 



(3) 



where r = (m c /m b ) 2 . The higher coefficient functions 6(r), c(r), ... are also perturbative series in a s . 
The expectation values of only two operators contribute to the \jm\ corrections: the Darwin term 
p% and the spin-orbital term p\ s . Due to the kinematic definition of the hadronic invariant mass 
M x , the general expression for the hadronic moments includes Mb explicitly: 



M n (M x ) = m 2 b n Yl 

1=0 



M B - m b 



' ' E nl (r) + a ni (r)— + b nl (r)^ + c nl (r)^ 

+d nl (rA + s nl (r)^ + ...}. (4) 

llli, Illy, 



It is possible to re-express the heavy quark masses, toq, in the above equations, in terms of the 
meson masses, M Hq , through the relation [6]: 

M HQ = m Q + A + 4^ + Pl + f LS 2 ~ Pll +0 (±) . (5) 
2m Q 4m 2 Q \ m Q ) 

The use of these expressions introduces an explicit dependence of the non-local correlators contribut- 
ing to Pni- In the notation of [7], p\ x corresponds to linear combinations of 7^_ 4 . 

Here, we employ the following two formalisms. The first one is based on the kinetic running 
masses, mg(/i), and non-perturbative parameters, introduced in [8]. No charm mass expansion is 
assumed. The second formalism employs quark pole masses and the flW and £>W meson mass 
relations. It represents a useful reference, as it has been already adopted in several studies. 

Contributions through O(a 2 /3 ) [9,10] and 0(\/m\) [11-14] to the moments are available. De- 
pending on the formulation adopted, the number of parameters involved at this order ranges from 
six to nine. Some of these parameters, like m b and A 2 ^ A*g/3, are relatively well known. Others, 
notably those which appear at 0(1/ ml), are virtually unknown. 

2.1 The mb(fj,), m c (/i) and Formalism 

The running kinetic quark masses m b (p) and m c (/j.) are considered here as two independent parame- 
ters. Apart from /J^(p) and defined as expectation values in the actual B meson, there are two 
l/m\ parameters, p 3 D and p\ s . The effect of p\ s turns out to be numerically small. In Eqs. (3) and 
(4) the mass ratio r is given by (m c (p) /m b (p)) 2 , and the b quark mass is understood as m b (p). The 
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Table 1: Numerical values of the coefficients in Eq.(3) evaluated at r=0.06 and m&(l GeV) =4.6 GeV. 
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Table 2: Numerical values of the coefficients in Eq.(4) evaluated at r=0.06 and m&(l GeV) = 4.6 GeV. 
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perturbative coefficients additionally depend on fi/rrib and the mass normalization scale \x is set at 
u, — 1 GeV. The functions in Eq. (3) are well-known parton expressions. The relevant coefficients 
are given in Table 1, for the central values of m b (l GeV) = 4.6 GeV and r ~ 0.06 obtained in our 
fit. Although we quote only the leading-order perturbative coefficients, we also include second-order 
BLM corrections in the analysis. Detailed expressions for the coefficients will be presented elsewhere. 

In the case of hadronic moments, we discard in Eq. (4) coefficients b n i, c n \ with I > 1, and d n i, 
s n i with l>0. The only non- vanishing E i0 coefficient is E 10 = r — M\,jm\. The value of the other 
coefficients, at r = 0.06, are listed in Table 2. Here we consider only 0(a s ) corrections and evaluate 
them using a s = 0.3. 



2.2 The A and Ai Formalism 

This widely used scheme results from the combination of the OPE with the HQET. Following the 
notation of Ref. [14], the moments are expressed in the following general form: 



Mr 



a + a\- 



a s (M £ 



at . A 



7T 



+ a 2 p ^ + h 



IT 



a s A ci Ai + c 2 A 2 + c 3 A 2 
+ b 2 — + 



7T M f 



Ml 



+ \d 1 \ 1 A + d 2 \ 2 A + d 3 A 3 + d 4 p 1 + d 5 p 2 + ^ d b+i % \ +o(^^-) 



(6) 



where k = n and k = 2n for leptonic and hadronic moments, respectively, while ao = for hadronic 

moments. Mb = 5.3135 GeV is the spin-averaged B meson mass. The second order BLM corrections 1 

are expressed in terms of (3q = 11 — 2/3n/, where we take rif = 3. The coefficients cii,bi,Ci,di are 

given in Table 3 for the first three leptonic, Mi^^(Ei) , and hadronic moments, Mi )2 ^(Mx)- In the 

1 The terms 0{a 2 s (3c)) and 0{a s K) are not available for the third hadronic moment. In our analysis we employ an 
estimate and the related uncertainty is included in the fit. 
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Table 3: Numerical values of the coefficients in Eq.(6). 
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leptonic case the a± coefficients agree with ref. [12], while the coefficients of the first two hadronic 
moments agree with [13,14]. Details of the derivation will be presented elsewhere. 

The non-perturbative parameters in Eq.(6) are related to those in Section 2.1 by the following 
relations, valid up to 0(a s ): 

Vir = -M ; t*a = 3A 2 + ; p D = pi ; p LS = 3p 2 • (7) 

777.5 TUb 

Perturbative corrections introduce a significant numerical difference between the parameters in the 
two schemes. At fx — 1 GeV: 

2 2 

A ~ M B - m b (l GeV) - \~ ^° - 0.26 GeV ; — Ai ~ pl(l GeV) - 0.17 GeV 2 . (8) 

A well known problem of this formalism is the instability of the perturbative series, due to the 
use of the pole quark masses. Large higher order corrections are however expected to cancel in the 
relation between physical observables, as long as all observables involved in the analysis are computed 
at the same order in a s [6, 15]. We also note that, as a consequence of the HQET mass relations for 
the mesons, the intrinsic expansion parameter in Eq.(6) is l/M^, rather than 1/Mb- The convergence 
of this expansion has been questioned, in view of indications [16, 17] that the matrix elements % of 
some non-local operators could be larger than that expected from dimensional estimates. 

3 Fits and Results 

This analysis is based on the preliminary Delphi measurements [3,4] of the first three moments of 
the hadronic mass and charged lepton energy, summarised in Table 4. Owing to the large boost of B 
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Table 4: Preliminary DELPHI results for the three leptonic and hadronic moments. 



Moment 


Result (stat) (syst) 


Mi (Ec) 
M 2 (E e ) 
M 3 (E e ) 


(1.383 ± 0.012 ±0.009) GeV 
(0.192 ± 0.005 ±0.008) GeV 2 
(-0.029 ± 0.005 ±0.006) GeV 3 


M 1 (M X ) 
M 2 (M X ) 
M 3 (M X ) 


(0.534 ± 0.041 ± 0.074) GeV 2 
(1.226 ± 0.158 ± 0.152) GeV 4 
(2.970 ± 0.673 ± 0.478) GeV 6 



hadrons in Z° — > bb events, the acceptance of the analyses can be extended down to the start of the 
lepton energy spectrum, making their theoretical interpretation more direct. The results correspond 
to B% and B~ mesons decays only. 

We perform a x 2 fit to these six moments, using the two theoretical frameworks discussed above. 
In the fit we also impose some additional constraints derived from independent determinations. 

In the kinetic mass scheme, we fit the full set of parameters: m b (l GeV), m c (l GeV), /i 2 , 
p 3 D and p 3 LS . We impose /4=(0.35±0.05) GeV 2 [16] and p| 5 =(-0.15±0.15) GeV 3 . Two mass 
constraints have also been applied: m b (l GeV) = (4.57±0.10) GeV [18], and, to be conservative, 
m c (l GeV) = (1.05±0.30) GeV. The most stringent is that on m b (l GeV). It must be noted that this 
constraint is largely equivalent to that derived from the first moment of the photon energy spectrum 
in b — > S7 in other studies [2]. Results are obtained for a s (m b ) = 0.22 ± 0.01 and are shown in 
Table 5. In order to study the effect of the bounds on m bjC introduced, the fit has been repeated 
unconstrained. Results are consistent, although the accuracy on the masses degrades. In particular 
we find mb(l GeV) = (4.61±0.15) GeV. It is interesting to observe that the mass constraints applied 
are of the scale of the fit sensitivity. Also, the central values of the heavy quark masses are in good 
agreement with independent determinations [18,19]. 

In the alternative approach based on pole masses, the fit extracts A, Ai, A2, pi and p 2 . We fix 
% = and impose two constraints from M B * — M B and M D * — M D which effectively reduce by two 
the number of free parameters. The results are given in Table 6. 

Projections of the constraints from the six moments in the m b -pl and m b -p 3 D planes are shown 
in Fig. 1 and those in the A-Ai and A-p x planes in Fig. 2. The xV n -d.f. of the fits is 0.96 and 0.35 
in the two formulations. Since the contributions proportional to p LS in the moment expressions are 
numerically suppressed, the fit is only marginally sensitive to its size and the result is determined by 
the constraint applied. By removing this, the fit would give p| s =(-1.0±0.7) GeV 3 . 

In contrast, the value of the leading 1/m 3 correction (parameterised by p 3 D or p\) can be deter- 
mined with satisfactory accuracy and its range agrees with theoretical expectations [16]. 

Systematic uncertainties due to ranges of residual parameters which have been fixed and missing 
terms in the expansions have been estimated. For the running mass formalism we propagate the 



5 



Table 5: Results of fit for the m b (fi), m c (u) and /i 2 (/i) formalism. 



Fit 


Fit 


Fit 


Syst. 




Parameter 


Values 


Uncertainty 


Uncertainty 




m b (l GeV) 


4.59 


± 0.08 


± 0.01 


GeV 


m c (l GeV) 


1.13 


± 0.13 


± 0.03 


GeV 


GeV) 


0.31 


± 0.07 


± 0.02 


GeV 2 


P?> 


0.05 


± 0.04 


± 0.01 


GeV 3 


Table 6: Results of fit for the A-Ai formalism. 


Fit 


Fit 


Fit 


Syst. 




Parameter 


Values 


Uncertainty 


Uncertainty 




A 


0.40 


± 0.10 


± 0.02 


GeV 


Ai 


-0.15 


± 0.07 


± 0.03 


GeV 2 


A 2 


0.12 


± 0.01 


± 0.01 


GeV 2 


Pi 


-0.01 


± 0.03 


± 0.03 


GeV 3 


P2 


0.03 


± 0.03 


± 0.01 


GeV 3 



uncertainty on a s and evaluate the effect of removing the BLM corrections from the lepton moments. 
In this scheme that is a small effect and higher order perturbative corrections are expected to be under 
control. Dimensional estimates suggest that \jm\ effects do not exceed the present experimental 
resolution. Other systematic uncertainties will be addressed in a dedicated publication. 

For the A-A x formalism we take the effect of % = (0.0±0.50) 3 GeV 3 , a s = 0.22±0.01 and we also 
estimate the effect of the missing corrections to third moments as M^(0. 001 ±0.0005) f3 (a s /rr) 2 and 
Mf (0.003 ± 0.003) A/ M B a s /n. 

The fit was also repeated using only the first two moments, leaving free m b (l GeV), /i 2 (l GeV) 
and A, Ai, respectively. The other parameters were fixed to the central values obtained in the full fit. 
Results agreed with those from the full fit. In particular, the values of A = 0.42±0.07(stat.) GeV and 
Ai = (— 0.17±0.05(stat.)) GeV 2 agree with the recent result reported by the Cleo Collaboration [1], 
which uses the first moments of the charged lepton energy to obtain A = 0.39 ± 0.07 GeV and 
Ai = (-0.25 ±0.05) GeV 2 . 

There are several facets of these results to be looked at. One interesting piece of information 
comes from the correlation between m c and m h extracted from the fit. It corresponds to m c (I GeV) = 
1.63 x (m b (l GeV) — 3.91). Therefore a competitive value of the charm mass can be extracted from 
a precise determination of m b . Using, for instance, rrib(l GeV) = (4.60 ± 0.05) GeV would give 
m c (l GeV) = (1.13 ± 0.09) GeV. This can be compared to the present typical lattice uncertainties 
which range between 50 and 120 MeV [20]. 

In the running mass scheme, the expansions of Eq.(5), for the B and D mesons are not used in 
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Figure 1: The projection of the constraints of the six measured moments on the rrib(l GeV)-/i^.(l GeV) 
(left) and m&(i GeV)-p% (right) planes. The bands correspond to the total measurement accuracy 
and are given by keeping all the other parameters at their central values. The ellipses represent the 
1 a contours. 



the fit. It is therefore possible to test a posteriori the consistency of the meson mass expansion by 
comparing the A values obtained in the two cases. We find A(B) — A(D) = — 0.086±0.092. This is 
also a test of the size of the non-local terms. 

In both approaches, the OPE predictions for the six moments, computed with the available 
precision, have a common intersection in the multi-parameter space and the quality of the fit is 
good. Within the present experimental accuracy, we therefore do not see the need to introduce 
higher order terms to establish agreement with the data. In particular, the first leptonic and hadronic 
moments are highly correlated and depend on nearly the same combination of heavy quark masses. 
Fixing this from Mi(Mx), one finds M\{E() = 1.377 GeV which agrees well with the measured 
value of (1.383±0.015) GeV. This provides a non-trivial consistency check of the OPE. The overall 
agreement represents both a test of the theory and suggests constraints on the size of the l/m^ terms 
and of other missing corrections. Similarly, the observed agreement strongly supports the validity of 
quark-hadron duality in the B decay shape variables. 

At present the achieved experimental resolution matches the available theoretical accuracy. With 
more precise data soon becoming available, it is important to improve the latter, particularly for 
higher hadronic moments. One way to improve the convergence of the heavy quark expansion could 
be to employ different kinematic variables. We propose to consider J\f x = M\ — 2AEx, where 
Mx and Ex are the hadronic mass and energy and A a fixed mass parameter. Choosing A near 
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Figure 2: The projection of the constraints of the six measured moments on the A-Ai (left) and A-pi 
(right) planes. The bands correspond to the total measurement accuracy and are given by keeping 
all the other parameters at their central values. The ellipses represent the 1 a contours. 



Mb — GeV) ~ 0.65 GeV, suppresses terms with I > 1 in Eq. (4) and results in a better 
convergence of higher moments [21]. The use of this variable should be feasible at B factories, where 
the kinematics allows an accurate reconstruction of both the mass and energy of the hadronic system 
in s.l. B decays. 

3.1 Implications for |V^| 

The value of |V c b| obtained from the total s.l. decay width depends on the OPE parameters extracted 
above. We discuss now the implications of our results for |V c b|, using the input parameters given 
in Table 7. The uncertainties on the BR(6 — > Xt~v) have been increased compared to ref. [22] for 
not using the heavy quark forward-backward asymmetries in the LEP global electroweak fit and to 
account for the ±15% uncertainty on the equality of s.l. partial width of b baryons and mesons. 

Table 7: Input values used for the determination of \V c b\. 



Measurement 


Value [22] 


6-hadron lifetime 
BR(6 -> Xe~v) 
BR(6 -> XjTv) 


(1.564 ±0.014) ps 
(10.59 ±0.31) % 
(0.17 ±0.05) % 
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The inclusive s.l. decay width has been calculated through second order in perturbative QCD. 

Second order BLM corrections were obtained in [23], all-order BLM terms are available from [24], 

whereas second-order non-BLM corrections have been estimated in [25] . Non-perturbative corrections 

start at order 0(l/mf) [11] and 0(l/m\) corrections have also been calculated [7]. Electroweak 

corrections have also been taken into account [26]. 

The determination of |V c b| and the contributions of the various parameters in the kinetic mass 
scheme is described in [27]. An approximate formula which displays the dependence on the different 
parameters is: 



|V cb | = |V cb | c 



- 0.65 (m 6 (l) - 4.6 GeV/c 2 ) + 0.40 (m c (l) - 1.15 GeV/c 2 ) 
+0.01 (y4 - 0.4 GeV 2 ) + 0.10 (pfj - 0.12 GeV 3 ) 
+0.06 (/4 - 0.35 GeV 2 ) - 0.01 (p| 5 + 0.17 GeV 3 ) 



(9) 



A detailed discussion of the theoretical uncertainties on jV^I goes beyond the scope of this pa- 
per. Here we focus on the uncertainty arising from the heavy quark masses and non-perturbative 
parameters determined in the fit. It is evaluated using the full fit error matrix which leads to ±1.5%. 
There is an additional uncertainty coming from the limited accuracy of the theoretical expressions 
which have been used. We take the range m^/2 < ji' < rrib for the scale ji' at which a s is evaluated 
and find a ±1% effect 2 . In summary, we obtain: 

|V cb | = 0.0419 x (l ± 0.016| meas ± 0.0151,., ± 0.010| peri ) , (10) 

where the first uncertainty reflects the accuracy on the s.l. width determination. 

The expression for the inclusive b s.l. width in the pole mass scheme is known to the same order. 
The fit results have been used to obtain: 

|Vcb| = 0.0413 x (l ± 0.016| meas ± 0.017| /it ± 0.006^ ± 0.021| pert ) . (11) 

Again, the first two uncertainties correspond to the decay width measurement and to the fitted 
parameters, respectively. The third uncertainty refers to the ^ =14 parameters which have been 
varied within the range (0 ± (0.5) 3 ) GeV 3 . The uncertainty from the truncation of the perturbative 
QCD series is again estimated by varying the scale at which a s is evaluated between m b /2 and Im^- 
Here the perturbative uncertainty is larger and reflects the slower convergence of the perturbative 
series when the pole mass scheme is employed. 



4 Conclusions 

The values of the heavy quark masses have been determined, together with the leading 1/m 2 and 

1/mf parameters, from a fit to the first three moments of the charged lepton energy and hadronic 

incorporating the third-order BLM correction suppresses this scale dependence. Combining Refs. [24] and [28], 
we find the third-order BLM correction to T s i(b^>c) to be s=a — 50(a s /7r) 3 in this scheme. This increases \V c b\ by 1% 
for n' — irib, and leaves it nearly unchanged, compared to two loops, for ^' = mb/2. 
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mass spectra in s.l. decays, measured in a preliminary analysis of the Delphi data. The absence 
of a charged lepton energy cut in the analysis makes the heavy quark expansion more reliable and 
allows us to include higher moments. We have adopted two different formalisms: one based on low- 
energy running quark masses, which does not rely on a l/m c expansion, and the other on pole quark 
masses. The constraints from the six moments agree well and the size of the dominant l/m\ term 
has been found to be compatible with theoretical estimates. The fit is largely insensitive to non-local 
correlators and to the spin-orbital operator. 

Propagating the ranges of the OPE parameters to the determination of | V c b\ reduces the theoretical 
uncertainty due to the 1/mf corrections below 2%. Furthermore, the use of a fit changes the nature of 
these uncertainties and partly removes the arbitrariness arising from estimates based on parameters 
ranges. 
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Note 

During the final stage of this work, a new analysis of s.l. moments has appeared [29]. There are 
several differences with our approach, but the results are qualitatively consistent with our findings. 
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